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ABSTRACT 

Context. Solar Shack-Hartmann wavefront sensors measure differential wavefront tilts as the relative shift between images from dif- 
ferent subapertures. There are several methods in use for measuring these shifts. 

Aims. We evaluate the inherent accuracy of the methods and the effects of various sources of error, such as noise, bias mismatch, and 
blurring. We investigate whether Z-tilts or G-tilts are measured. 

Methods. We test the algorithms on two kinds of artificial data sets, one corresponding to images with known shifts and one corre- 
sponding to seeing with different r$. 

Results. Our results show that the best methods for shift measurements are based on the square difference function and the absolute 
difference function squared, with subpixel accuracy accomplished by use of two-dimensional quadratic interpolation. These methods 
measure Z-tilts rather than G-tilts. 

Key words, atmospheric effects - instrumentation: adaptive optics - site testing 



1. Introduction 

High-resolution observational astronomy with telescopes oper- 
ated on the ground relies on methods for combating the effects 
from turbulence in the Earth's atmosphere. Turbulence varies 
with time scales on the order of ms, mixing air that varies in 
temperature and therefore in refractive index. The wavefronts of 
the light from astronomical objects, which are flat when they 
enter the atmosphere, are distorted by the optically active atmo- 
sphere, which causes image motion and blurring in the images 
collected through the telescope. Measurements of the wavefronts 
allow us to correct for these effects, both in real time and in post 
processing. 

There are several different methods for measuring in real 
time the variation of the wavefront phase of the light over the 
telescope pupil. For solar telescopes, the one exclusively used is 
the Shack-Hartmann (SH; |Shack & Platt|1971| > wavefront sen- 
sor (WFS). This is true for the adaptive optics (AO) systems of 
present 1 -meter class, high-resolution solar t elescopes (|Rim mele 
& Radick|1998l|Scharmer et al.|2000[|2003l |von der Liihe et al 



2003; Rimmele et al. 2004). SH wide-field WFSs are used for 



characterizing the atmospheric turbulence above existing and fu- 
ture solar telescopes ( Waldmann 2007 ; Waldm ann et al.| ]2008 
IScharmer & van Wer khoven 2 01 0[ >. 

An SH WFS estimates the wavefront from measurements of 
the local wavefront tilt in several regions of the pupil called sub- 
pupils. The local tilts manifest themselves as image motion in 
the focal plane associated with each subpupil. For night-time 
telescopes, this image motion is (comparably) easily measured 
by tracking the peak of the image of a natural or artificial point 
source. As no point sources are available on the sun, one instead 
has to resort to measuring relative image motion by matching 
solar features in the different subpupil images. This can be done 



by finding the minimum of the degree of mismatch of the images 
as a function of image shift. This procedure requires solving two 
distinct subproblems: 1) computing the mismatch function as a 
function of image shift in steps of whole pixels and 2) finding 
the minimum of the mismatch function with subpixel precision. 
The latter involves interpolation between the whole-pixel grid 
points. This procedure for solar correlation tracking was appar- 



ently pioneered by Smithson & Tarbell ( |1977[ ) and | von der Ltihe 
( fT953] l. 

The SH WFSs developed for the AO systems currently in 
use at different solar telescopes use a variety of shift measure- 
ment methods, so do the sensors used for turbulence characteri- 
zation. This suggests that the choices of methods were based on 
personal preferences and technological momentum from other 
applications rather than on a thorough evaluation of the perfor- 
mance of these methods on relevant data. 

Several algorithms for measuring the mismatch and perform- 



ing the interpolation were investigated by Yi & Molowny Horas 
( T992J. However, because the purpose of their work was differ- 
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ent from ours (image remapping in post-processing), they used 
much larger fields of view (FOVs) than is currently practical for 
real-time solar SH WFS. 

As part of a project on characterization of the atmospheric 
turbulence using SH WFS, Waldmann] ( |2007| l and [Waldmann 
et al. (2008) address these problems for setups more similar to 
ours. The recent MSc thesis of Johansson] (2010) has expanded 
on Waldmann's work in this area. She implemented a number of 
different methods for image-shift measurements and tested them 
on synthetic data relevant to solar SH WFSs. Comparisons with 
their methods and results will appear throughout this paper. 

In this paper we investigate, by use of artificial data rele- 
vant to present solar SH WFSs, a number of different algorithms 
for measuring whole-pixel image shifts and interpolating to get 
subpixel accuracy. The algorithms are defined in Sect. [2] We in- 
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Table 1: Correlation algorithms (CAs). 



Acronym Name 



Mismatch, c,j, for shift i, j 



SDF 


Square Difference Function 




g Kt {x + i,y + j)f 


CFI 


Covariance Function, image domain 




■ g^(x + i,y + j) 


CFF 


Covariance Function, Fourier domain 




g(x t y)} ■ g*{w 2 • glS(x,y)}}(i, j) 


ADF 


Absolute Difference Function 




ftefC* + i,y + j)\ 


ADF 2 


Absolute Difference Function, Squared 


x,y 


-gnt(x + i,y + J)\) 



Notes, g denotes the Fourier transform and an asterisk as a superscript denotes complex conjugation, wi is a 2D Hamming window. The notation 
g is used for a version of g, where a fitted plane has been subtracted. 



vestigate the inherent accuracy of the algorithms in Sect. [3] by 
testing them on identical images with known shifts, as well as 
the influence of noise and variations in intensity level. In Sect. [4] 
we use the best methods on images formed through artificial see- 
ing and evaluate the performance for different seeing conditions. 
We discuss our results in Sect. 



2. Algorithms 

2.1. Correlation algorithms 

In Table[T]we define five different correlation algorithms (CAs), 
which we use to calculate the image mismatch on a grid of inte- 
ger pixel shifts i, j. The names and acronyms of the CAs are also 
given in the table. These mismatch values make a matrix c with 
elements c,j. A coarse estimate of the image shift, (Sx, <5y), is 
then given by the indices corresponding to the grid position with 
the minimum mismatch value, (i m i n , y m i n ). This shift should be 
sought within a maximum range in order to reduce the number 
of false matches to other parts of the granulation pattern. The 
are then used to refine this estimate to 



algorithms in Sect. 2.2 
subpixel accuracy. 

Perhaps most straight-forward is the SDF algorithm, in 
which one calculates the mismatch in a Least Squares (LS) 
sense. 

Subtracting the intensity mean, expanding the square in the 
SDF algorithm, and retaining only the cross term gives twice the 
covariance (with negative sign), which is the basis of the follow- 
ing two algorithms. The CFI algorithm calculates the covariance 
in the image domain. It is the one being used for the Dunn Solar 
Telescope system (Rimmele & Radick 1998). The correlation 
coefficient differs from the covariance only by division with the 
standard deviations of the two images, so methods based on the 
former (e.g., Waldmann et al.|2008|l should give results similar 



to those of CFI. 

The covariance can also be calculated in the Fourier domain, 
taking advantage of the Fast Fourier Transform (FFT). The CFF 
algorithm was developed by von der Luhe| ( T983) for an image 
stabilization system and is used today in the KAOS AO imple- 
mentation used at the Vacuum Tower Telescope ( |von der Liihe| 
et al. 2003| ). For small images, such as those involved in SH 
WFS calculations, one can expect errors from wrap-around ef- 
fects. That is, because of the assumption of periodicity implicit 
in finite-size Fourier transforms, for large shifts, structures in one 
image are not matched with structures at the shifted location but 
with structures shifted in from the opposite side of the image. 



|Knutsson et al.| ((2005) derived a method based on the 
Fourier spectrum of the correlation function (correlation spec- 
trum phase). It is supposed to give results similar to the CFF 
method but with some accuracy sacrificed for speed. We have 
not investigated this method. 

The CFF method requires apodization of the images, i.e., the 
multiplication of a window function that reduces ringing effects 
from the discontinuities caused by the Fourier wrap-around. 
When not explicitly stated otherwise, we use a Hamming win- 
dow written as 

W2(x,y) = wi(x) ■ wi(y), (1) 



where w\ is the 1 -dimensional Hamming window (Enochson & 
|Otnes|1968| >, 



w\{x) — a + (a 



(2) 



See 



where a = 0.53836 and N is the linear size of the window, 
also Sect. l3.3.5l below. 

The ADF algorithm is fast because it can be calculated very 
efficiently in CPU instructions available for many architectures, 

1 p003l > use ADF for the IR 



particularly for 8-bit data. Keller et al. 



AO system of the McMath-Pierce so 



ar telescope. So do IMiura] 
et al.| ( |2009[ ) in their recently presented AO system used for the 
domeless solar telescope at Hida Observatory. 

Because the shape of the ADF minimum does not match the 
assumption of a parabolic shape implicit in the subpixel algo- 
rithms, squaring the ADF correlation values leads to an improve- 
ment. This adds a completely negligible computational cost to 
that of ADF, as squaring is done after summing. In fact, be- 
cause it does not move (Zmin, jmm)> only the at most nine grid 
points used for subpixel interpolation (see Sect. 2.2 1 have to be 
squared. This ADF 2 method was developed by G. Scharmer in 
1993 for use in the correlation tracker systems of the former 
Swedish Vacuum Solar Telescope (Shand et al. |1995[ l and is in 
use in the AO and tracker systems of the Swedish 1 -meter Solar 
Tel escope (SST; |Scharmer et al.|2 000 2003}. 

|Smithson & TarbellT"( fl977| l showed that a linear trend in 
intensity shifts the covariance peak from the correct position. 
Therefore one should subtract a fitted plane from both g and g re f 
before applying the CAs. However, the granulation data used 
in our simulations have negligible trends and for the difference 
based algorithms (ADF, ADF 2 , and SDF), a consistent bias in the 
intensity level cancels automatically. We therefore saved com- 
puting time in our tests by limiting the pre-processing of the data 
to only subtracting the mean values, and only when calculating 
CFI and CFF. 
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Table 2: Interpolation algorithms (I A) or subpixel minimum finding algorithms. 



Acronym 


Name 


Location of minimum 












r„ _ (/ „ A • — / c , A 1 /? 












fir, — (7 e, A ■ — 9/ e„ A . 4- / c 


j>/)/2 


2LS 


2D Least Squares 


\xi » yi ) 


where 


- (7? , V - ,W2 












as = «S(,i>i - 2(Si,o)i + 


>,)/2 










«6 = - 5-1,1 - Sl,-1 + 5- 


■i,-i)/4 










'«2 = 0i,o - «-i,o)/2 












«3 = (Ji,o - 2^o,o + *-i,o)/2 




2QI 


2D Quadratic Interpolation 


(X2,y2) 


where • 


«4 = («o,i - *o,-i)/2 












«5 = (*o,i - 2%o + so,-i)/2 












«6 = (Sl,l - 5-1,1 - «1,-1 + «- 


-1.-0/4 


1LS 


ID Least Squares 




using a 2 , 03, «4, a s from 2LS above 




1QI 


ID Quadratic Interpolation 


(x u yi) 


using a 2 , 03, «4, a 5 from 2QI above 





Notes. Notation in the table refers to the elements, Sy, of the submatrix s of the correlation matrix centered on its sample minimum, ((minOmin), 
the expression for a second order polynomial of two variables, f(x, y) = ay + a 2 x + a^x 2 + a 4 y + a 5 y 2 + a&xy, and the location of its ID minimum, 
x i — ! min ~~ a 2 /2a^; y\ = j m i„ — 124/205 and its 2D minimum, x 2 = z m j„ + {2a 2 as — a 4 a^)j{a 1 b — 4a^as)', y 2 = jmm + (2fl3fl4 — a 2 a(,)/(a^ — 40305). The 
notation (•>, denotes averaging over index i e {-1,0, 1}. 



2.2. Subpixel interpolation 

The methods in Sect. |2.1| are responsible for making a coarsely 
sampled 2D correlation function with a reasonable shape. The 
interpolation algorithms (IAs) in this section then have to find 
the minimum with better accuracy than given by the sampling 

grid^ 

November & Simon] (p"988 ) found that interpolation meth- 



ods should be based on polynomials of degree 2. Methods based 
on polynomials of higher degree systematically underestimate 
the shift for small displacements, while first degree polynomials 
give a systematic overestimation. The algorithms we consider 
can all be described as fitting a conic section, 

f(x,y) — a\ + CI2X + a^x 2 + a 4 y + asy 2 + a$xy, (3) 

to the 3x3-element submatrix s of c centered on the sample min- 
imum of c with elements 



,,J+Jm 



I, J 



-1,0,1. 



(4) 



The interpolated shift vector, (6x, 5y), is the position of the min- 
imum of the fitted function, (x m i n ,y m [ n ). The algorithms differ 
in the number of points used and whether the fitting is done in 
2D or in each dimension separately. The definitions, names, and 
acronyms of the different methods are given in Table [2] 

The 1QI method is based on numerical ID derivatives using 
the center row (column) of s ( November| 1 986 1 . It is equivalent 
to a LS estimate using only the center row (column) of s. It does 
not use the corner elements of s. 

The 1LS algorithm consists of, separately for the x and y 
directions, fitting a ID polynomial of degree 2 to all the ele- 
ments in s. This is equivalent to the procedure of Waldmann 



(2007), projection of the data onto the axes (i.e., summing the 
rows (columns)), and doing LS fitting on the result. Waldmann 
uses Lagrange interpolation but this is mathematically equiva- 
lent to a LS fit. 

The 2QI algorithm was derived by |Yi & Molowny Horas 



( 1992 their Eq. (9)) as an extension of the 1QI algorithm. 

The 2LS algorithm is based on expressions for the conic sec- 
tion coefficients derived by Waldmann (2007 ). Johan sson] ( |2010| l 



found that Waldmann's expression for one of the 2LS coeffi- 
cients is missing a factor 2. We are using the corrected expres- 
sion. 

|Waldmann| ( |2007| l compared four different subpixel meth- 
ods: 1LS, 2LS, and a 2D fit to a Gaussian. He found that the 
Gaussian fit gave the best results and the polynomial method 
worked almost as well. Because the Gaussian fit is more compu- 
tationally heavy, he adopted the latter. Johansson ( |2010 i tested 
the Gaussian fit and got results comparable to Waldmann's but 
much better results for the polynomial fits. Based on Johansson's 
results, we have not evaluated Gaussian fits as a method for sub- 
pixel interpolation. 

Miura et al. ( |2009 1 use a method for subpixel interpolation, 
where they find the centroid of a spot generated by inverting and 
then clipping the mismatch measured with ADF. This method is 
not considered here. 



3. Algorithm accuracy 

The goal of this experiment with artificial data is to establish the 
accuracy of the algorithms themselves, for granulation images 
(GI), which are identical except for a known shift and detector 
imperfections such as noise and bias. In reality, the images will 
be different because the atmospheric turbulence not only shifts 
the images but also smears them differently. Those effects will 
be addressed in Sect. |4]below. 

3.1. Artificial data 

For this experiment, the images should be identical except for 
a known shift, defined to subpixel precision without introduc- 
ing errors that are caused by the re-sampling needed for sub- 
pixel image shifting. We therefore start with a high-resolution 
image, degrade the resolution, shift it by a known number of 
whole high-resolution pixels, and then down-sample it to the SH 
image scale. 

Specifically, we use a 2000x2000-pixel, high-quality SST 



la 



G-band image with an image scale of O'.'041/pixel, see Fig 
The image was recorded on 25 May 2003 by Mats Carlsson et 
al. from ITA in Oslo and corrected for atmospheric turbulence 
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(a) (b) 

Fig. 1: Granulation image used for simulation experiments, (a) 430.5 nm (G-band) SST diffraction limited granulation image (GI). 
2000x2000 pixels, image scale (X'041/pixel. (b) Degraded to resolution of a 9.8 cm subpupil, shifted by a whole number of pixels, 
and re-sampled to 200x200 0741-pixels. 



effects by use of Multi Frame Blind Deconvolution (Lofdahl 
|2002) ). We degraded it to 9.8 cm hexagonal (edge to edge) sub- 
pupil resolution at 500 nm. This degraded image was shifted by 
integer steps from to 20 times the high-resolution pixels, as 
well as in steps of 10 pixels from 30 to 60. The so degraded and 
shifted GIs were box-car compressed by a factor 10 to 200x200 
pixels of size 0741. This procedure gives images with known 
subpixel shifts, 5x and 6y, without any re-sampling, except for 
the compression. Figure lb shows a sample compressed image. 

The data with 0-20 high-resolution pixel shifts were made 
to test subpixel accuracy, while the 30-60-pixel shifts are for 
testing linearity with larger shifts. 

The diffraction limited resolution, /L/D m \, * 1" at 500 nm, 
corresponds to > 2 pixels. This means subpixel accuracy corre- 
sponds to super-resolution accuracy. 

The resulting images had more contrast than the real data 
from our SH WFS, so some bias was added to change the RMS 
contrast to ~3% of the mean intensity. The resulting images were 
stored in two versions, with and without Gaussian noise with a 
standard deviation of O.S'^Q The digitization noise of a 12-bit 
camera is insignificant compared to the Gaussian noise but may 
be significant for an 8-bit camera. We do not include the effects 
of digitization in our simulations. 



3.2. Processing 

The 200x200-pixel FOV is much larger than the FOV of the 
SH WFS, which allows the use of many different subfields in 
order to get better statistics. Centered on each of 17x17 grid 
positions, subimages, g, of size 16x16 or 24x24 pixels, were 
defined. The subimages defined in the unshifted reference image, 
g re f, were larger in order to accommodate a shift range limited to 
+8 pixels along each axis direction, except for CFF, which uses 
two images of equal size. Note also that for CFF, the size of the 
correlation matrix is limited by the subimage size. For 16x16- 



1 This noise level corresponds to photon noise from a CCD with 40 
ke" full well. For a WFS set up to be running all day, the noise level 
would be larger when the sun is at low elevation and the image therefore 
darker. However, it is likely that the performance is then limited by other 
effects than noise, such as image warping from anisoplanatism. 



pixel subfields, this limits the range to +6 pixels (in reality to 
even less). 

The different sizes of g, 16x16 and 24x24 pixels, have two 
purposes: 1) We want to see how a change in size affects some 
of the methods and 2) we will compare CFF using 24x24 pixels 
with the other methods using 16x16 pixels. If the image geome- 
try on the detector accommodates a 24x24-pixel g K f, then it can 
also accommodate 24x24-pixel g subfields if no oversize refer- 
ence image is needed. 

We measured the shifts with each combination of CAs in 
Sect. 12. H and IAs in Sect. 12.21 We do this with and without noise 
and with and without multiplying the reference image by 1.01, 
giving an approximate 1% bias mismatch. The bias mismatch 
sensitivity is investigated because it is known to be a problem 
with the ADF and ADF 2 methods. 

For each real shift, 5x lea i (and 5y lea i, we will simplify the 
notation by referring to both axis directions with x when possi- 
ble), we calculate the mean and standard deviation, <x, of the 
measured shifts, <5;t measure d, by use of the outlier-robust statis- 
tics based on Tukey's Biweight as implemented in the IDL 
Astronomy User's Library (Landsm an|1993 i. 

After removing 4<x outliers, we fit the data to the relationship 



^measured = P0 + Pl&Xreal + Pi SHl(a (5x re al), 



(5) 



where a = 2n/ 1 pixel, by use of the robust nonlinear LS curve 
fitting procedure MPFIT ( |More|[T978} |Markwardt||2009"| l. The 
linear coefficient, p\, of these fits is listed in the tables below. 

3.3. Results 

In Fig. [2] the SDF, ADF 2 , CFI, and CFF CAs are compared for 
the case of no noise and no bias mismatch, using the 2QI IA and 
24x24-pixel subfields. The errors are a mix of systematic and 
random errors. The SDF and ADF 2 algorithms give a tighter cor- 
relation between the true and measured image shifts than the CFI 
and CFF algorithms. November & Simon (1988]) found similar 
undulating effect of small errors near whole and half pixel shifts 
calculated by CFF and larger errors in between. Balle steros et al.| 
(1996) make a similar observation with ADF. 

The ADF and CFF methods produce many outliers, which 
is the main reason we need the robust statistics mentioned in 
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Fig. 2: Errors in measured shifts vs. real shifts for four different CFs, zero noise N, and zero bias mismatch B. Mean values and +lcr 
error bars calculated with robust statistics. The dashed line in the CFF panel corresponds to pi — 1, where p\ is fitted to Eq. d5J. 



Sect. |3.2| A very small fraction of these outliers, on the order 
10~ 5 of all cases for CFF, zero for all other CAs, are caused 
by the correlation matrix minimum falling on an outer row or 
column. The rest are caused by false minima in the correlation 
matrix, which happen to be deeper than the real minimum. This 
can happen occasionally if a secondary minimum is located on 
or near a grid point, while the real minimum is between grid 
points. The effect is more severe for ADF because of the more 
pointy shape of its minimum. For CFF, minima corresponding to 
large shifts are attenuated by two effects: apodization lowers the 
intensity away from the center of the subimage and the digital 
Fourier transform wraps in mismatching structure from the other 
side of the subimage. This can lead to detection of false minima 
corresponding to small shifts. 

The dashed line in the CFF panel represents the linear part of 
the fit to Eq. PJ; its slope is p\ — 1. The CFF method systemati- 
cally underestimates large shifts. Comparison between the fitted 
line and the mean values shows a slight nonlinearity in the CFF 
method, making the underestimation worse for larger shifts. 

The deviation of p\ from unity and the undulations are sys- 
tematic errors. The former can be calibrated while the latter mix 
with the random errors represented by the error bars. 

Figure [^illustrates the effects of noise and of bias mismatch 
for SDF. Here we limit the plots to 6r = (&e 2 eal + ^ eal ) 1/2 < 
2 pixels. The blue line corresponds to the fit of Eq. The am- 
plitudes of the undulations, pi, is ~0.01 pixel. Noise does not 
appear to change the undulating errors significantly but it does 



make the random errors more dominant. Bias mismatch has a 
similar effect and overwhelms the differences due to noise. It 
also appears to change a minute systematic overestimation of 
the shifts to an equally small underestimation. 



3.3.1. Tables 

The complete results of these simulations are presented in 
Tables HI] 

While the standard deviations in the plots are for each real 
shift individually, and therefore do not include the undulations, 
the standard deviations in the tables are for intervals of real shifts 
and therefore do include the undulations. The measured shifts 
might be expected to have near-Gaussian distributions, which 
is true for ADF 2 and SDF and mostly for CFI. ADF often has 
complicated, multi-peak distributions. All CFF distributions are 
double-peaked and/or asymmetrical. 

We give the results for all shifts, 6r, but also separately 
for small shifts (< 1 pixels and < 2 pixels, resp.), medium 
shifts (shifts between 3 and 5 pixels in length), and large shifts 
(> 5 pixels). The small shifts are relevant to AO performance in 
closed loop. The large shift results tell us something about per- 
formance in open loop, which is relevant for wavefront sensor 
calibration, site testing, image restoration, and when trying to 
close an AO loop. 
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Fig. 3: Results using The SDF/2QI methods with 24x24-pixel subfields. Small shifts, 5r S5 2 pixels, and different combinations of 
noise, N, and bias mismatch, B, as indicated in the labels of each tile. Error bars are +lcr. The blue line is a fit to Eq. (|5J. 



3.3.2. Identical images 

In Table [3] we show the performance of the different methods 
with noise free data. 

We begin by noting that in many cases the errors for large 
shifts are smaller than those for small shifts. This may seem 
counterintuitive but it is simply a consequence of using only 
whole pixel shifts for the larger shifts. The IAs have smaller sys- 
tematic errors on the grid points, see the undulations in Fig. [2] 

The ADF CA produces many 4<x outliers. It is not surprising 
that it gives much worse results than ADF 2 , because the two CAs 
by definition share the location of the whole pixel minimum but 
the ADF minimum does not have the parabolic shape assumed 
by the IAs. 

SDF and ADF 2 are clearly better than the CFI and CFF meth- 
ods. CFF can compete for the very smallest shifts (< 1 pixel) if 
we use 24x24 pixels for CFF and 16x16 pixels for the other CAs 
(which is the most fair comparison, since the optics need to ac- 
commodate 24x24 pixels in order to get the oversize reference 
image for the non-CFF CAs). 

Except for CFF, the errors appear more or less independent 
of the magnitude of the shifts. CFF deteriorates significantly in 
three ways at larger shifts: The number of outliers increases, the 
random error increases, and the shifts are systematically under- 
estimated, as shown by the nonunity slopes. The former two ef- 
fects can be explained by the assumption of periodicity of the 
digital Fourier transform. The signal for large shifts is diluted by 



mismatching granulation shifted in from the opposite end of the 
FOV. 

For SDF and ADF 2 , the 2D IAs are clearly better than the ID 
ones. 2QI is marginally better than 2LS. The best results with 
SDF and ADF 2 are an error RMS of less than 0.02 pixels, cor- 
responding to 07008. Increasing the subfield size from 16 to 24 
pixels squared reduces the error by approximately 30%. 



3.3.3. Noisy images 

Adding 0.5% noise to the images give the results shown in 
Table [4] The errors grow but the behavior is similar to the zero- 
noise case. 

The best results with SDF and ADF 2 16x16 are an error 
RMS of less than 0.03 pixels, corresponding to 07012. CFF with 
24x24 gives similar results for small shifts. As in the case of zero 
noise, increasing the subfield size from 16 to 24 pixels squared 
reduces the error by approximately 30%. 

With 24x24 pixels, SDF and ADF 2 give results similar to the 
16xl6-pixel zero-noise case. Surprisingly, the number of out- 
liers for ADF is significantly reduced by the added noise. But cr 
increases more than for ADF, so this may be an effect of making 
the error distribution more Gaussian. 
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Fig. 4: Window functions, 24x24 pixels. Left: Hamming win- 
dow. Right: Flat-top window. 



3.3.4. Bias mismatch 

Subtracting a fitted plane (or just the mean intensity as in our 
experiment, see Sect. 



2.1 



from each subimage removes mis- 
matches in intensity bias for CFI and CFF We did not do this 
for the difference based CAs (SDF, ADF, and ADF 2 ), where 
a consistent bias cancels. However, if there is a bias mismatch 
between the images this cancellation is not effective. Such bias 
mismatch could come from, e.g., variations in a thin cloud layer 
in the case when g re f is not from the same exposure as g, or small 
drifts in the pupil location on the SH, causing variations in the 
light level from the outermost subpupils. 

We re-processed all the data after multiplying the reference 
image by 1.01, introducing a 1% bias mismatch with the other 
images. In Tables [5] and [6] we show the results for the differ- 
ence based CAs. As expected, the CFI and CFF results did not 
change from the ones in Tables [3] and [4] and they are therefore 
not repeated. 

For small shifts, the SDF and ADF 2 results are now worse 
than for CFF, even when comparing the methods using the same 
image size. SDF is more robust against bias mismatch than 
ADF 2 . There is now no real difference between ADF and ADF 2 , 
possibly indicating that the parabolic shape of the ADF 2 correla- 
tion function is destroyed by the bias mismatch. The SDF error 
RMS is 0.7 pixels for 16x16 and 0.4 pixels for 24x24, which 
makes it the best method for large shifts. 

With a bias mismatch, 0.5% noise added to the images does 
not significantly change the results for the difference based CAs. 
We conclude that bias mismatch should be removed in pre- 
processing by subtraction of the intensity mean. 



3.3.5. Window function 

The CFF method requires apodization, i.e., the multiplication 
of the subimage (after subtracting the fitted plane) by a win- 
dow function. The intention is to reduce ringing effects from the 
discontinuities caused by Fourier wrap-around. For the results 
above, we used a 2D Hamming window. In Tables [7J|8] we show 
the results when we instead use a window with the center ~50% 
of the area flat, and a taper only in the pixels outside this area. 
This is similar to the window used by Waldmann (2007). See 
Fig.|4]for the two types of window functions. 

The results are mixed. The errors are larger but there are 
fewer outliers. The results are similar with and without noise. 
The systematic underestimation is reduced, resulting in slopes 
closer to unity. Based on these data we cannot say which win- 
dow is better, they would have to be evaluated specifically for 
any new situation where one wants to use CFF. The important 
result for our purposes is that the comparison between SDF and 
ADF 2 vs CFF does not depend on the window function used. 



4. Image shift as a measure of wavefront tilt 

For wavefront sensing, the quantities that really need to be mea- 
sured are the average wavefront gradients at the positions cor- 
responding to the subapertures. One assumes that a shift in im- 
age position corresponds perfectly to the average gradient of the 
wavefront across the subaperture. However, in addition to local 
gradients, continuous wavefront aberrations across the telescope 
aperture also result in local wavefront curvature. Is the assump- 
tion valid anyway and how good is it for different seeing condi- 
tions, as quantified with Fried's parameter ro? 

Using Kolmogorov statistics for different ro without any as- 
sumption of partial correction by an AO system makes the statis- 
tics from this experiment relevant to open-loop WFS. I.e., sys- 
tems for measuring seeing statistics (e.g., Scharmer & van 



|Werkhoven|20 1 Q\ and the capture phase for AO systems, but not 
necessarily AO systems in closed loop. 



4.1. Artificial data 

The setup corresponds to a filled 98 -cm pupil (like the SST) with 
85 subapertures, each 9.55 cm edge-to-edge. In Fig. 5a we show 
one sample Kolmogorov wavefront phase. The superimposed 
pattern shows the geometry of 85 hexagonal microlenses cir- 



cumscribed by a telescope pupil. Figure 5b demonstrates the lo 



cal plane approximation implicit in SH wavefront sensing, while 
Fig. 5c shows just the tilts. 

The exact geometry does not matter for the results reported 
here, since our tests do not involve the step where wavefronts are 
reconstructed from the shift measurements. However, the geom- 
etry discussed is along the lines planned for the next generation 
SST AO system, and in that respect motivates the particular sub- 
aperture size and shape investigated. 

We generated 100 wavefront phases, (pi, following 
Kolmogorov statistics. For making each simulated phase 
screen, the following procedure was used. 1003 random num- 
bers were drawn from a standard normal distribution, scaled 
with the square root of the atmospheric variances and used as 
coefficients for atmospheric Karhunen-Loeve (KL) functions 
2-1004. We used KL functions based directly on the theory of 
Fried ( 1978 1, as implemented by |Dai| ( 1995 1. These modes are 
numbered in order of decreasing atmospheric variance, and the 
exact range of indices is motivated by KL] being piston and 
KL 1( )05 a circular mode, starting a higher radial order. Figure 5a 
shows a sample wavefront masked with the pattern of the 85 
hexagonal microlens geometry. 

For all simulations, we used a wavelength of 500 nm and 
a telescope aperture diameter of D te \ = 98 cm. In order to 
cover a range of different seeing conditions, we scaled these 
wavefront phases to different values of Fried's parameter, ro e 
{5,7, 10, 15,20} cm, by multiplying with (D tei /r ) 5/6 . Using the 
same wavefront shapes scaled differently like this, the perfor- 
mance for different values of ro should be directly comparable. 

For each random wavefront, separately scaled to each value 
of ro, and for each subpupil defined by a microlens, we generated 
an image by convolving the GI in Fig. la with a PSF based on 
the subpupil and the local wavefront phase. 

We want to examine the effect of using different subfield 
sizes. Increased subfield size can be used in different ways: 
Either one can change the image scale, so the same amount of 
granulation fits in the FOV but in better pixel resolution. Or one 
can keep the original image scale so more granulation fits in the 
FOV. (Or something in between.) Therefore we make images 
at three different image scales by box-car compressing them by 
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(a) 



(b) 



(c) 



Fig. 5: Sample wavefront and the subpupil pattern of 85 hexagonal microlenses. The circle represents the 98 cm SST pupil and the 
gaps between adjacent microlenses have been exaggerated for clarity. The geometry is actually such that there is no space between 
the microlenses. (a) The 1003-mode simulated wavefront phase, (b) Piecewise plane approximation, (c) Local tilts only. 



three different integer factors, 7, 10, and 13. The resulting image 
scales are 0729/pixel, 0741/pixel, and 0753/pixel. 

To summarize: the images we have generated were down- 
graded to the resolution of the subpupil, shifted by the local 
wavefront tilt and also somewhat blurred by the local wavefront 
curvature, the latter in particular for data with small ro. A bias 
was then added to make the RMS contrast of the granulation 
pattern approximately 3% of the mean intensity. 

4.2. Processing 

For each subfield size, noise, and ro, relative shifts and tilts 
were calculated for 490,000 randomly selected pairs of subpupil 
images. We operate on image pairs corresponding to subpupils 
from different random wavefronts, so results are not influenced 
systematically by spatial correlations. 

For the shift measurements, we use all the CF methods from 
Sect. [2] except ADF, applied to the two images in a pair. We use 
only the 2QI method for subpixel interpolation. 

Shift measurements were calculated twice, with and without 
1% noise added to the images. We increased the noise level from 
the 0.5% used in Sect. !3.3.3l to make the effect of noise clearer 
and thus allow better comparison of different methods. We do 
not investigate bias mismatch in this experiment. The conclusion 
from Sect. [3] is that bias mismatch should be compensated for 
before applying the shift measurement methods. 

For comparison with the shift measurements, <5x s hift, we fit- 
ted Zernike tip and tilt to each wavefront, within each hexagonal 
subpupil. The relative wavefront tilt for an image pair is the dif- 
ference between the Zernike tilts for the two subpupils in the 
pair. These relative tip/tilt coefficients in radians, a x , are con- 
verted to image shift, 

Sxftx = 2 ^ a x (6) 
7irD tel 

where r is the image scale in rad/pixel. We calculate the robust 
statistics of the resulting shifts, remove Act outliers, and fit the 
data to the relationship 

foshift = Po + Pi Sxm + P2 sin(a 5x t iit), (7) 
where a = 2n/l pixel. Compare Eq. Q. 

4.3. Results 

Table [9] shows the results from images with image scale 
0741/pixel and different subfield sizes, NxN pixels with N = 16, 



24, and 36. In seconds of arc, this corresponds to 6'.'56, 9784, 
and 14776 squared. The results using the additional image scales 
0729/pixel and 0753/pixel and a single array size, N = 24, are in 
Table [TO] 

In addition to the tabulated results, we calculated the follow- 
ing: The standard deviation of dxnu varies linearly with r 5 ^ 6 as 
expected from Kolmogorov statistics: cr ti i, = 2.68, 2.02, 1.50, 
1.07, 0.84 pixels = 1710, 0783, 0762, 0744, 0/34 for r = 5, 7, 
10, 15, 20 cm, resp. Correlation coefficients calculated after re- 
moval of outliers are very high: they round to 1 .00 in all cases, 
except CFF and CFI at N = 16 for which they are 0.98-0.99. 



4.3.1. Failures 

If the errors were normal distributed, using a Act limit for defin- 
ing outliers should give a failure rate of 0.0063%. With almost 
5 • 10 5 samples we expect the normal distribution to be well re- 
alized but the actual failure rates are larger. All CAs show an 
excess but by far it is the CFF results that suffer from the highest 
number of outliers, particularly for small subfields and small ro. 
CFI is also slightly worse than SDF and ADF 2 . 

All failure rates are ^ 10%. They are > 2% for ro = 
5 cm with all methods when using the smallest image scale, 
0729/pixel, and for small FOVs and r = 5, 7 cm (and CFI 
ro = 5 cm, noise, N = 16). Figure [6] illustrates how failure rates 
sometimes increase with noise, sometimes decrease. With noise, 
the failure rate decreases with large FOV in arcsec (exception: 
N = 36 and ro = 7 cm). The variation is less systematic without 
noise. 



4.3.2. Shift errors 

We now plot some of the results in pixels rather than seconds of 
arc so we can compare with the results in Sect. [3] 

Figure|7 shows the easiest case: the largest subfields, 36x36 
pixels, and the best seeing, ro = 20 cm, resulting in the smallest 
shifts. The ADF 2 results are so similar to the SDF results that we 
only show the latter of those two methods. 

The errors in Fig.|7]are consistent with the ones in Fig. [2] We 
recognize the undulations but we also see a linear trend, not only 
for CFF. SDF and CFF: cr en . x 0.025 matches undulations plus 
error bars (CFF because for 20 cm we are dominated by small 
shifts). CFI: cr en * 0.045 is actually slightly better than in Fig. [2] 
These results (as opposed to the ones in Fig. [8} are comparable 
because 20 cm is much larger than the subpupil size so the im- 
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ages are not blurred by local phase curvature. We cannot directly 
compare the noisy data, because here we have 1 % noise, in the 
earlier experiment we used 0.5%. 

Linear fit offsets are small, \po\ < 0.002 pixels for all CAs, 
except CFI that has \po\ < 0.016 pixels. The undulations fit sines 
with an amplitude of p2 ~ 0.01 pixels in almost all cases, con- 
sistent with Figs. [2] and [3] where we used identical images with 
known shifts. The exceptions are cases when both N and ro are 



small. In the latter cases the fits to the sine function sometimes 
give as small amplitudes as po ~ 0.002 pixels or less. 

For shifts smaller than 0.2 pixels, the overall overestimation 
turns into underestimation. For CFF, the slope and the undula- 
tions work in the same direction, making the underestimation 
larger for the smaller shifts. 

Figure [8] shows the most demanding test: the smallest sub- 
fields, 16x16 pixels, and the worst seeing, ro = 5 cm, resulting 
in the largest shifts and the images most blurred by local phase 
curvature. Note much larger dispersion and error in p\ for CFF. 

4.3.3. Z-tilts and G-tilts 

The wavefront tilt measured over a subpupil is by necessity an 
approximation because in reality the tilts vary over the sub- 
pupil. In the night-time literature, two kinds of tilts are dis- 
cussed (e.g., Tokovinin 2002 ). G-tilts correspond to averaging 
the wavefront Gradient, which is mathematically equivalent to 
measuring the center of Gravity of the PSF However, because 
of noise and asymmetrical PSFs this can never be realized in 
practice. Windowing and thresholding the PSF gives measure- 
ments that are more related to Z-tilts, corresponding to Zernike 
tip/tilt and the location of the PSF peak. When interpreting the 
measurements, we need to know what kind of tilts are measured 
by our methods. 

The simulated tilts are implemented as coefficients to the 
Zernike tip and tilt polynomials. So if the shifts measure Z-tilts, 
the expected p\ is by definition 



E[pi | Z-tilt] EE 1. 



(8) 



In order to derive the expected p\ for G-tilts, we use formulas 
given by |Tokovinin| ( |200"2) . The variance of the differential im- 
age motion can be written as 



KA 2 r: 5/3 D- 1 ' 3 



(9) 



where D is the subpupil diameter and K is a number that depends 
on the kind of tilt. The expected p\ should be equal to the ratio 
of <x d for G-tilts and Z-tilts, i.e., 



E[pi | G-tilt] = tJk g /K z * 0.966, 



(10) 



where we used K = K G = 0.340 for G-tilts and K - K z - 0.364 
for Z-tilts, which is asymptotically correct for large separations 
between subpupils. This corresponds to a 3.4% difference in tilt 
measurements with Z-tilts giving the larger numbers, regardless 
of wavelength, seeing conditions, and aperture size. For smaller 
separations, E[p\ | G-tilt] is even smaller and depends somewhat 
on whether the shifts are longitudinal or transversal (parallel or 
orthogonal to the line separating the apertures). 

In our results for all the CAs except CFF, 1.007 ^ p\ $ 
1.010. This means they overestimate Z-tilts systematically by 
^1% and G-tilts by ^ 4.6%. This result appears to be robust 
with respect to noise, subfield size, image scale, and seeing con- 
ditions. 

The CFF p\ values confirm the underestimation of the im- 



age shift found in Sect. 3.2 They depend mostly on the size of 
the subfield but also to some extent on ro and noise. CFF p\ im- 
proves with larger FOV in arcsec, regardless of number of pixels 
and ro. For small FOVs, CFF p\ is smaller than the G-tilts slope 
of 0.966, grows with FOV size to cross G-tilt at ~14". For the 
other CAs, p\ is independent of FOV size. CFF slope improved 
with subfield size. Is it the size in pixels or the size in arcsec that 
is important? Figure|9]makes it clear that the size in arcsec is the 
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important parameter and that with even larger FOVs we should the important parameter is the number of 1" resolution elements 
expect to measure something closer to Z-tilts. We conjecture that that fit within the FOV. 
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4.3.4. RMS errors 

The RMS errors, cr elT , are calculated after removing the outliers 
and after subtracting the linear fit. In spite of having more out- 
liers removed from the calculations, the errors are worse for CFF 
than for the other methods. 

Compared to the noise-free data, adding 1% image noise 
significantly increases <x err . In many cases it also decreases the 
number of fails (all cases with N = 24, 36, most cases with 
N = 16). This probably means the noise makes the error dis- 
tribution more Gaussian. 

<x en decreases with increasing ro, indicating that the algo- 
rithms perform better in good seeing, as expected. But are the 
results worse in bad seeing because of local phase curvature (i.e., 
smearing) or just because the shifts are larger? The relative mea- 
sure, cr err /cr t iit, does not decrease as much with r ( > for zero-noise 
data and actually tends to increase for noisy data. So the latter 
should not be the major effect. 

How do the errors depend on the FOV size in pixels and in 
seconds of arc? Figure 10 visualizes the following results from 
the tables. Without noise, the CFI and CFF results improve with 
N as well as the FOV size in arcsec with constant N (particularly 
from 6'.'9 to 978). SDF (and ADF 2 , the results are so similar that 
we show only SDF) errors are more or less independent of FOV 
size, whether in pixels or in arcsec. With 1% noise, however, for 
all methods, the errors improve only slightly or grow when N 
is constant and we increase the FOV size in arcsec (particularly 
from 7" to 13"). When N is increased from 16 to 24 with almost 
constant 7" FOV, there is a significant improvement. There is 
also significant improvement when N is increased from 24 to 36 
(although FOV in arcsec is increased too). So for the RMS error, 
FOV in pixels is more important than FOV in arc seconds with 
noisy data. 
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Blue: ro = 7 cm. N x A^-pixel subfields. Diamonds (O): im- 
age scale 0741/pixel, varying N = {16,24,36} as labeled in the 
plots; Down triangles (v): 0729/pixel, N = 24; Up triangles (A): 
0753/pixel, N = 24. Symbols below the dotted line correspond to 
data without noise, above the line with 1 % noise. Note different 
vertical scales. 



5. Discussion and conclusions 

We have evaluated five different correlation algorithms and four 
different interpolation algorithms in two experiments with arti- 
ficial data. Among these are the algorithms in use for the AO 
systems installed at the major high-resolution solar telescopes 
today. The experiment in Sect. [3] examined the inherent perfor- 
mance of the methods with identical images and known shifts, 
and introduced also the effects of noise and a mismatch in in- 
tensity level (bias) between images. We used image contrast and 
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noise levels that resemble the setup at the SST. In Sect. [4] we 
introduced also local wavefront curvature and different seeing 
conditions. 

Based on the results of both experiments, we recommend 
SDF and ADF 2 (in that order) for calculating the correlation 
functions. They give significantly smaller random errors and 
more predictable systematic errors than the competing methods. 
The SDF results are marginally better so if computational speed 
is not an issue, use SDF. But ADF 2 may, by virtue of speed, 
be preferable. For subpixel interpolation, 2QI and 2LS perform 
better than the one-dimensional interpolation algorithms. 

It is clear from Sect. [3] that bias mismatch has a strong and 
negative effect on the performance of the difference based meth- 
ods. We have not established how much bias mismatch can be 
tolerated. Our conclusion is that it should be compensated for 
in pre-processing, before the shift measurement methods are ap- 
plied. The simulated bias mismatch comes from multiplication 
with 1.01 but is compensated for by subtraction for CFI and 
CFF Based on the fact that the latter methods give identical re- 
sults with and without the bias mismatch (to the number of digits 
shown in the table), we conclude that it does not matter whether 
the mismatch is removed by multiplication or by subtraction. 

|Waldmann| ( |2007| ) considers only a single correlation algo- 
rithm (CFF) but tries a few different interpolation algorithms: 
1LS, 2LS, and a Gaussian fit. He finds a Gaussian fit to give 
the better results but the almost as good 1LS is less complicated 
so he uses that. On the other hand, Waldmann et al. ( 2008 1 try 



CFF, CFI (almost: correlation instead of covariance) and ADS, 
but use only 1LS for interpolation. They find CFF best and ADF 



worst. Trying all combinations of those methods, Johansson 



( |2010| l confirmed that Gaussian and 1LS perform similarly to- 
gether with CFF. She also found a similar performance for SDF 
and ADF 2 with 1LS and Gaussian, but significantly better per- 
formance when interpolated with 2QI and 2LS. Based on this, 
we did not try a Gaussian fit here because of its greater compu- 
tational cost. 

We have demonstrated that with the recommended methods, 
we can measure shifts in identical, noise-free 24x24-pixel im- 
ages with an RMS error of 0.012 pixels, corresponding to 07008 
in our setup. With realistic image noise RMS of 0.5% the errors 
increase by 40% to 0.017 pixels. Increasing the noise RMS to 
1 % gives unacceptably high measurement errors in bad seeing, 
particularly with small subimages. Increasing the subfield size 
from 16 to 24 pixels squared reduced the error by approximately 
30%. 

One method is not necessarily the best choice both for open 
loop and for closed loop. In open loop, performance for large 
shifts and predictable systematic errors are important. However, 
CFF can compete with SDF and ADF 2 in closed loop, where 
wavefronts are already compensated for and the CFF random er- 
rors are small. The underestimation of the shifts can be reduced 
by using a larger FOV and can also be compensated for with a 
different servo gain. Based on our results, it is difficult to find a 
reason to ever use the CFI method. 

We have found that SDF, ADF 2 , and CFI measure Z-tilts 
with a 1% systematic overestimation, rather than G-tilts. For the 
FOV sizes relevant to solar SH WFS, the CFF method severely 
underestimates the tilts but it is likely that for larger FOVs, also 
CFF measures Z-tilts. 

Johansson ( 2010| l processed also images together with a ref- 



erence image that was slightly distorted geometrically (com- 
pressed image scale in one direction and expanded in the other) 
and found that it could significantly affect the results. We have 
not examined this effect here. A relevant further test involving 



anisoplanatic effects would require geometrical distortions and 
differential blurring whose shape and magnitude can be calcu- 
lated from realistic atmospheric turbulence profiles. 

We end by emphasizing that small details in processing may 
have large effects on the results. Examples from our evaluation 
include ADF versus ADF 2 and the choice of apodization for 
CFF. As in so many other situations, it is important to keep track 
of implementation details and other tricks, so they are not lost 
when upgrading SH software. 
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Table 3: Simulation results, no bias mismatch, no noise. 



CA 


IA 




All 


Sr < 1 


Sr< 2 


3 <S 


r<5 


Sr > 5 


■ Pi 


it ( ni 


Fail (%} 


cr (pix) 


Fail (%) 


cr (pix) 


Fail (%) 




Fail (%) 


cr (pix) 


Fail (%) 


SDF 


1LS 


0.037 


1.0 


0.042 


1.5 


0.041 


1.4 


0.028 


0.3 


0.027 


0.3 




SDF 


1QI 


0.033 


1.8 


0.039 


2.5 


0.038 


2.4 


0.023 


0.6 


0.022 


0.6 




SDF 


2LS 


0.020 


0.3 


0.021 


0.3 


0.020 


0.3 


0.020 


0.1 


0.020 


0.3 




SDF 


2QI 


0.017 


0.3 


0.017 


0.2 


0.017 


0.3 


0.016 


0.2 


0.016 


0.3 




CFI 


1LS 


0.114 


0.8 


0.119 


0.5 


0.118 


0.7 


0.109 


0.9 


0.106 


0.9 




CFI 


1QI 


0.108 


0.7 


0.112 


0.5 


0.112 


0.6 


0.103 


0.9 


0.100 


0.8 




CFI 


2LS 


0.109 


0.7 


0.109 


0.4 


0.110 


0.6 


0.109 


1.0 


0.105 


0.9 




CFI 


2QI 


0.101 


0.6 


0.102 


0.4 


0.102 


0.5 


0.101 


0.9 


0.099 


0.9 




CFF 


1LS 


0.112 


5.0 


0.060 




0.078 




0.187 


3.0 


0.295 


30.1 


0.86 


CFF 


1QI 


0.111 


5.0 


0.063 




0.079 




0.180 


3.0 


0.280 


30.1 


0.86 


CFF 


2LS 


0.091 


5.1 


0.032 




0.056 




0.167 


3.2 


0.269 


30.8 


0.86 


CFF 


2QI 


0.085 


5.1 


0.031 




0.054 




0.156 


3.2 


0.251 


30.9 


0.86 




1 T Q 




Q Q 
y.y 


0.049 


14.3 


0.049 


13.7 




l.y 


0.037 


2.6 




ADF 


101 


0.042 


33.0 


0.046 


39.3 


0.047 


39.1 


0.031 


22.4 


0.027 


20.6 




ADF 


2LS 


0.030 


23.1 


0.033 


25.1 


0.033 


25.1 


0.022 


19.2 


0.022 


17.5 




ADF 


201 


0.028 


42.8 


0.032 


49.4 


0.032 


49.7 


0.017 


30.1 


0.016 


27.6 




ADF 2 


1LS 


0.038 


0.9 


0.044 


1.3 


0.044 


1.3 


0.029 


0.3 


0.028 


0.3 




ADF 2 


lOI 


0.035 


1.5 


0.041 


2.2 


0.041 


2.1 


0.025 


0.5 


0.023 


0.5 




ADF 2 


2LS 


0.022 


0.3 


0.023 


0.4 


0.023 


0.4 


0.021 


0.1 


0.021 


0.1 




ADF 2 


20T 


0.019 


0.1 


0.019 


0.1 


0.019 


0.1 


0.018 




0.018 






SDF 


1LS 


0.025 


1.2 


0.029 


1.7 


0.029 


1.6 


0.019 


0.4 


0.019 


0.5 




SDF 


1QI 


0.024 


1.8 


0.028 


2.6 


0.028 


2.5 


0.018 


0.6 


0.017 


0.7 




SDF 


2LS 


0.014 


0.7 


0.014 


0.9 


0.014 


0.8 


0.013 


0.3 


0.013 


0.2 




SDF 


2QI 


0.012 


0.8 


0.012 


0.9 


0.012 


0.9 


0.012 


0.5 


0.011 


0.4 




CFI 


1LS 


0.070 


0.3 


0.073 


0.2 


0.072 


0.2 


0.067 


0.3 


0.064 


0.4 




CFI 


1QI 


0.066 


0.2 


0.070 


0.2 


0.069 


0.2 


0.062 


0.3 


0.060 


0.4 




CFI 


2LS 


0.065 


0.3 


0.067 


0.1 


0.066 


0.2 


0.065 


0.5 


0.062 


0.6 




CFI 


2QI 


0.061 


0.2 


0.062 




0.062 


0.1 


0.061 


0.5 


0.058 


0.5 




CFF 


1LS 


0.053 


0.8 


0.042 




0.047 




0.055 


0.5 


0.083 


4.5 


0.94 


CFF 


1QI 


0.053 


0.8 


0.043 




0.048 




0.052 


0.5 


0.079 


4.4 


0.94 


CFF 


2LS 


0.040 


0.6 


0.017 




0.025 




0.062 


0.4 


0.092 


3.6 


0.94 


CFF 


2QI 


0.036 


0.7 


0.017 




0.024 




0.055 


0.5 


0.082 


3.8 


0.94 


\ 1 ' 1 


1 T ^ 


PI P19Q 




0.032 


41.3 


0.033 


41.7 


O 09 1 


99 7 


0.019 


20.8 




ADF 


1QI 


0.027 


49.6 


0.030 


58.7 


0.031 


59.4 


0.015 


31.6 


0.012 


29.0 




ADF 


2LS 


0.016 


47.9 


0.019 


54.3 


0.019 


55.0 


0.009 


33.8 


0.010 


31.1 




ADF 


2QI 


0.014 


57.3 


0.018 


66.9 


0.019 


67.6 


0.009 


37.0 


0.008 


34.0 




ADF 2 


1LS 


0.026 


1.2 


0.029 


1.8 


0.029 


1.6 


0.020 


0.4 


0.019 


0.5 




ADF 2 


1QI 


0.025 


1.7 


0.029 


2.5 


0.029 


2.3 


0.018 


0.6 


0.017 


0.7 




ADF 2 


2LS 


0.015 


0.6 


0.015 


0.8 


0.015 


0.7 


0.014 


0.1 


0.014 


0.3 




ADF 2 


2QI 


0.013 


0.2 


0.014 


0.4 


0.014 


0.3 


0.013 


0.1 


0.012 


0.1 





Notes. See Tables[T|and|2]for the CA and IA acronyms. The "Fail (%)" columns show the percentage of 4cr outliers when not rounded to 0.0. The 
fitted slope, pi, is shown only when not rounded to 1.00. The FOV is 16x16 pixels in the upper part and 24x24 pixels in the bottom part. 
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Table 4: Simulation results, no bias mismatch, 0.5% noise. 



CA 


IA 




All 


Sr < 1 


Sr< 2 


3 <S 


r<5 


Sr > 5 


■ Pi 


it ( ni 


Fail (%} 


cr (pix) 


Fail (%) 


cr (pix) 


Fail (%) 




Fail (%) 


cr (pix) 


Fail (%) 


SDF 


1LS 


0.041 


0.6 


0.045 


1.0 


0.045 


0.9 


0.033 


0.2 


0.033 


0.2 




SDF 


1QI 


0.038 


1.0 


0.043 


1.4 


0.043 


1.3 


0.030 


0.3 


0.029 


0.2 




SDF 


2LS 


0.026 


0.1 


0.026 


0.1 


0.026 


0.1 


0.026 


0.1 


0.026 


0.1 




SDF 


2QI 


0.024 


0.1 


0.024 


0.1 


0.024 


0.1 


0.023 


0.1 


0.023 


0.1 




CFI 


1LS 


0.117 


0.7 


0.122 


0.5 


0.121 


0.7 


0.111 


0.9 


0.107 


0.8 




CFI 


1QI 


0.110 


0.7 


0.115 


0.4 


0.114 


0.6 


0.105 


0.8 


0.101 


0.8 




CFI 


2LS 


0.111 


0.7 


0.113 


0.4 


0.113 


0.6 


0.111 


0.9 


0.107 


0.9 




CFI 


2QI 


0.104 


0.6 


0.105 


0.4 


0.105 


0.5 


0.103 


0.9 


0.100 


0.8 




CFF 


1LS 


0.118 


5.0 


0.069 




0.085 




0.189 


3.1 


0.297 


30.0 


0.86 


CFF 


1QI 


0.117 


5.0 


0.072 




0.086 




0.182 


3.1 


0.281 


30.0 


0.86 


CFF 


2LS 


0.097 


5.1 


0.043 




0.064 




0.168 


3.3 


0.271 


30.7 


0.86 


CFF 


2QI 


0.092 


5.2 


0.043 




0.062 




0.158 


3.3 


0.252 


30.8 


0.86 




1 T Q 


ft CiAQ 


A ft 


0.053 


6.2 


0.053 


5.8 


n nAn 

U.U^fU 


n 7 
u. / 


0.038 


0.7 




ADF 


101 


0.055 


11.2 


0.059 


16.1 


0.059 


15.6 


0.047 


3.3 


0.044 


3.1 




ADF 


2LS 


0.039 


5.8 


0.041 


6.5 


0.041 


6.5 


0.033 


4.6 


0.032 


4.3 




ADF 


201 


0.041 


20.2 


0.044 


23.5 


0.045 


23.6 


0.033 


14.2 


0.032 


13.1 




ADF 2 


1LS 


0.043 


0.6 


0.048 


0.9 


0.048 


0.8 


0.035 


0.1 


0.034 


0.1 




ADF 2 


lOI 


0.041 


0.8 


0.046 


1.2 


0.046 


1.1 


0.032 


0.2 


0.031 


0.2 




ADF 2 


2LS 


0.028 


0.2 


0.029 


0.2 


0.029 


0.2 


0.027 


0.1 


0.028 


0.1 




ADF 2 


20T 


0.026 


0.1 


0.027 


0.1 


0.027 


0.1 


0.025 




0.025 






SDF 


1LS 


0.028 


0.7 


0.032 


1.0 


0.032 


0.9 


0.022 


0.2 


0.022 


0.3 




SDF 


1QI 


0.027 


0.9 


0.031 


1.3 


0.031 


1.3 


0.021 


0.3 


0.021 


0.4 




SDF 


2LS 


0.018 


0.1 


0.018 


0.2 


0.018 


0.2 


0.017 


0.1 


0.017 


0.1 




SDF 


2QI 


0.017 


0.1 


0.017 


0.1 


0.017 


0.1 


0.016 




0.016 






CFI 


1LS 


0.071 


0.3 


0.074 


0.2 


0.073 


0.2 


0.067 


0.4 


0.065 


0.4 




CFI 


1QI 


0.067 


0.2 


0.071 


0.2 


0.070 


0.2 


0.063 


0.3 


0.061 


0.3 




CFI 


2LS 


0.066 


0.3 


0.068 


0.1 


0.067 


0.2 


0.066 


0.5 


0.064 


0.5 




CFI 


2QI 


0.061 


0.2 


0.063 




0.063 


0.1 


0.061 


0.5 


0.059 


0.5 




CFF 


1LS 


0.057 


0.8 


0.047 




0.051 




0.059 


0.5 


0.086 


4.3 


0.94 


CFF 


1QI 


0.057 


0.7 


0.049 




0.053 




0.057 


0.5 


0.083 


4.1 


0.94 


CFF 


2LS 


0.045 


0.6 


0.026 




0.032 




0.065 


0.4 


0.093 


3.5 


0.94 


CFF 


2QI 


0.042 


0.6 


0.026 




0.031 




0.058 


0.4 


0.083 


3.6 


0.94 


\ 1 ' 1 


1 T ^ 




Ij.O 


0.040 


18.0 


0.041 


18.0 


n rn 1 


D.U 


0.029 


5.8 




ADF 


1QI 


0.038 


29.9 


0.042 


35.4 


0.042 


36.0 


0.029 


19.0 


0.027 


17.5 




ADF 


2LS 


0.027 


22.5 


0.029 


23.6 


0.029 


24.7 


0.021 


18.2 


0.020 


16.7 




ADF 


2QI 


0.026 


39.2 


0.029 


44.1 


0.030 


45.4 


0.018 


27.9 


0.017 


25.5 




ADF 2 


1LS 


0.029 


0.7 


0.033 


1.0 


0.033 


0.9 


0.023 


0.2 


0.023 


0.3 




ADF 2 


1QI 


0.029 


0.8 


0.032 


1.2 


0.032 


1.1 


0.022 


0.2 


0.022 


0.3 




ADF 2 


2LS 


0.019 


0.1 


0.020 


0.2 


0.020 


0.2 


0.018 




0.018 


0.1 




ADF 2 


2QI 


0.018 




0.018 




0.018 




0.017 




0.017 







Notes. See Tables[T|and|2]for the CA and IA acronyms. The "Fail (%)" columns show the percentage of 4cr outliers when not rounded to 0.0. The 
fitted slope, pi, is shown only when not rounded to 1.00. The FOV is 16x16 pixels in the upper part and 24x24 pixels in the bottom part. 
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M.G. Lofdahl: Evaluation of image-shift measurement algorithms for solar SH WFS 
Table 5: Simulation results, 1% bias mismatch, no noise. 



CA 


IA 




All 


Sr< 1 


Sr < 2 


3 < Sr < 5 


Sr> 5 




Fail (%\ 


tr (pix) 


Fail (%) 


a (pix) 


Fail (%) 


a- (pix) 


Fail (%) 


a (pix) 


Fail (%) 


SDF 


1LS 


0.087 


0.1 


0.092 


0.1 


0.090 


0.1 


0.082 




0.083 




SDF 


1QI 


0.080 


0.1 


0.084 


0.2 


0.083 


0.2 


0.074 




0.075 




SDF 


2LS 


0.079 


0.1 


0.081 


0.1 


0.079 


0.1 


0.079 


0.1 


0.081 




SDF 


2QI 


0.071 


0.1 


0.072 


0.1 


0.071 


0.2 


0.071 


0.1 


0.072 




a r*r-? 
\ i j t 


1 T <J 


n 1 ha 


w.z 


0.112 


0.2 


0.109 


0.2 


0.095 


0.2 


0.097 


0.1 


ADF 


101 


0.107 


0.3 


0.115 


0.3 


0.112 


0.3 


0.098 


0.3 


0.099 


0.3 


ADF 


2LS 


0.100 


0.2 


0.104 


0.2 


0.102 


0.2 


0.095 


0.2 


0.097 


0.1 


ADF 


2QI 


0.100 


0.3 


0.105 


0.2 


0.103 


0.3 


0.096 


0.4 


0.098 


0.4 


ADF 2 


1LS 


0.108 


0.1 


0.113 


0.1 


0.111 


0.1 


0.102 




0.104 




ADF 2 


lOI 


0.108 


0.1 


0.111 


0.1 


0.109 


0.1 


0.106 


0.1 


0.108 




ADF 2 


2LS 


0.100 


0.1 


0.103 


0.1 


0.100 


0.2 


0.100 


0.1 


0.102 




ADF 2 


2QI 


0.101 


0.2 


0.101 


0.2 


0.099 


0.2 


0.103 


0.2 


0.106 


0.1 


SDF 


1LS 


0.055 


0.1 


0.058 


0.1 


0.058 


0.1 


0.050 




0.049 




SDF 


1QI 


0.050 


0.1 


0.053 


0.1 


0.053 


0.1 


0.045 




0.044 




SDF 


2LS 


0.049 




0.049 




0.050 




0.047 




0.047 




SDF 


2QI 


0.043 




0.044 




0.044 




0.042 




0.041 




\ 1 1 \ 


1 T *C 


u.uo / 


yi.j 


0.071 


0.3 


0.071 


0.4 


0.059 


0.1 


0.058 


0.1 


ADF 


1QI 


0.069 


0.3 


0.074 


0.4 


0.074 


0.3 


0.062 


0.1 


0.060 


0.1 


ADF 


2LS 


0.063 


0.2 


0.065 


0.1 


0.065 


0.2 


0.059 


0.1 


0.058 


0.1 


ADF 


2QI 


0.064 


0.2 


0.066 


0.3 


0.066 


0.3 


0.060 


0.2 


0.058 


0.2 


ADF 2 


1LS 


0.068 




0.070 




0.070 




0.063 




0.062 




ADF 2 


1QI 


0.069 




0.070 




0.070 




0.066 




0.065 




ADF 2 


2LS 


0.062 




0.063 




0.063 




0.061 




0.060 




ADF 2 


2QI 


0.064 




0.064 




0.064 




0.064 




0.063 





Notes. See Tables[T]and|2]for the CA and IA acronyms. The "Fail (%)" columns show the percentage of 4cr outliers when not rounded to 0.0. The 
fitted slope, pi, is shown only when not rounded to 1.00. The FOV is 16x16 pixels in the upper part and 24x24 pixels in the bottom part. 



Table 6: Simulation results, 1% bias mismatch, 0.5% noise. 



CA 


IA 




All 


Sr< 1 


Sr < 2 


3 < Sr < 5 


Sr> 5 


a (pix) 


Fail (%) 


a- (pix) 


Fail (%) 


cr (pix) 


Fail (%) 


or (pix) 


Fail (%) 


a (pix) 


Fail (%) 


SDF 


1LS 


0.089 


0.1 


0.094 


0.1 


0.092 


0.1 


0.084 




0.086 




SDF 


1QI 


0.082 


0.1 


0.086 


0.2 


0.085 


0.1 


0.076 




0.077 




SDF 


2LS 


0.081 


0.1 


0.083 


0.1 


0.081 


0.1 


0.081 




0.083 




SDF 


2QI 


0.073 


0.1 


0.074 


0.1 


0.073 


0.1 


0.073 




0.074 




ADF 


1LS 


0.103 


0.2 


0.110 


0.2 


0.108 


0.2 


0.094 


0.1 


0.096 




ADF 


1QI 


0.103 


0.3 


0.111 


0.3 


0.109 


0.3 


0.093 


0.2 


0.094 


0.2 


ADF 


2LS 


0.098 


0.2 


0.103 


0.2 


0.101 


0.2 


0.094 


0.1 


0.096 


0.1 


ADF 


2QI 


0.096 


0.3 


0.101 


0.2 


0.099 


0.2 


0.091 


0.3 


0.092 


0.3 


ADF 2 


1LS 


0.107 


0.1 


0.112 


0.1 


0.110 


0.1 


0.101 




0.103 




ADF 2 


1QI 


0.103 


0.1 


0.107 


0.2 


0.105 


0.1 


0.100 


0.1 


0.102 




ADF 2 


2LS 


0.099 


0.1 


0.101 


0.1 


0.099 


0.1 


0.098 




0.101 




ADF 2 


2QI 


0.096 


0.1 


0.097 


0.2 


0.095 


0.2 


0.097 


0.1 


0.099 


0.1 


SDF 


1LS 


0.057 


0.1 


0.060 




0.060 


0.1 


0.052 




0.051 




SDF 


1QI 


0.052 


0.1 


0.055 


0.1 


0.055 


0.1 


0.047 




0.047 




SDF 


2LS 


0.051 




0.052 




0.052 




0.049 




0.049 




SDF 


2QI 


0.046 




0.046 




0.046 




0.044 




0.044 




ADF 


1LS 


0.067 


0.3 


0.071 


0.4 


0.072 


0.4 


0.060 




0.059 


0.1 


ADF 


1QI 


0.068 


0.3 


0.073 


0.5 


0.073 


0.4 


0.059 


0.1 


0.058 


0.1 


ADF 


2LS 


0.064 


0.2 


0.066 


0.2 


0.066 


0.2 


0.059 


0.1 


0.058 


0.1 


ADF 


2QI 


0.063 


0.3 


0.065 


0.3 


0.066 


0.3 


0.058 


0.2 


0.056 


0.2 


ADF 2 


1LS 


0.068 




0.071 




0.071 




0.063 




0.062 




ADF 2 


1QI 


0.066 




0.068 




0.068 




0.063 




0.061 




ADF 2 


2LS 


0.063 




0.063 




0.064 




0.061 




0.060 




ADF 2 


2QI 


0.061 




0.061 




0.062 




0.060 




0.059 





Notes. See Tables[TJand|2]for the CA and IA acronyms. The "Fail (%)" columns show the percentage of 4cr outliers when not rounded to 0.0. The 
fitted slope, pi, is shown only when not rounded to 1.00. The FOV is 16x16 pixels in the upper part and 24x24 pixels in the bottom part. 
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M.G. Lofdahl: Evaluation of image-shift measurement algorithms for solar SH WFS 
Table 7: Simulation results, no bias mismatch, no noise. Flat top window for FFT. 



CA 


IA 




All 


6r < 1 


Sr < 2 


3 < 6r < 5 


Sr>5 


1 i 


cr (pix) 


Fail (%) 


tr (pix) Fail (%) 


a (pix) Fail (%) 


c (pix) 


ran y /c) 


cr (pix) 


rail y/c) 


CFF 


1LS 


0.146 


2.0 


0.076 


0.105 


0.193 


1.6 


0.305 


11.6 


0.95 


CFF 


1QI 


0.141 


1.9 


0.076 


0.103 


0.183 


1.5 


0.293 


11.6 


0.95 


CFF 


2LS 


0.136 


2.0 


0.063 


0.095 


0.187 


1.7 


0.295 


11.8 


0.95 


CFF 


2QI 


0.131 


2.0 


0.062 


0.092 


0.176 


1.5 


0.282 


11.7 


0.95 


CFF 


1LS 


0.077 


0.4 


0.044 


0.063 


0.082 


0.5 


0.122 


1.7 


0.97 


CFF 


1QI 


0.074 


0.3 


0.044 


0.061 


0.077 


0.5 


0.116 


1.6 


0.97 


CFF 


2LS 


0.072 


0.3 


0.034 


0.054 


0.082 


0.4 


0.123 


1.5 


0.97 


CFF 


2QI 


0.068 


0.3 


0.032 


0.052 


0.076 


0.3 


0.115 


1.4 


0.97 



Notes. See Tables[T|and|2]for the CA and IA acronyms. The "Fail (%)" columns show the percentage of 4<r outliers when not rounded to 0.0. The 
fitted slope, pi, is shown only when not rounded to 1.00. The FOV is 16x16 pixels in the upper part and 24x24 pixels in the bottom part. 
Table 8: Simulation results, no bias mismatch, 0.5% noise. Flat top window for FFT. 



CA 


IA 




All 


Sr< 1 


Sr<2 


3 < 6r < 5 


Sr> 5 


Pi 


cr (pix) 


Fail (%) 


cr (pix) Fail (%) 


cr (pix) Fail (%) 


cr (pix) 


Fail (%) 


cr (pix) 


Fail (%) 


CFF 


1LS 


0.148 


2.0 


0.079 


0.107 


0.194 


1.6 


0.307 


11.8 


0.95 


CFF 


1QI 


0.144 


2.0 


0.079 


0.105 


0.185 


1.5 


0.295 


11.7 


0.95 


CFF 


2LS 


0.139 


2.0 


0.066 


0.098 


0.188 


1.6 


0.295 


12.0 


0.95 


CFF 


2QI 


0.133 


2.0 


0.065 


0.095 


0.177 


1.4 


0.283 


11.9 


0.95 


CFF 


1LS 


0.078 


0.3 


0.046 


0.064 


0.083 


0.5 


0.123 


1.6 


0.97 


CFF 


1QI 


0.076 


0.3 


0.046 


0.063 


0.078 


0.4 


0.117 


1.5 


0.97 


CFF 


2LS 


0.073 


0.3 


0.036 


0.056 


0.084 


0.3 


0.124 


1.3 


0.97 


CFF 


2QI 


0.069 


0.3 


0.035 


0.054 


0.078 


0.3 


0.116 


1.2 


0.97 



Notes. See Tables[TJand|2]for the CA and IA acronyms. The "Fail (%)" columns show the percentage of 4cr outliers when not rounded to 0.0. The 
fitted slope, pi, is shown only when not rounded to 1.00. The FOV is 16x16 pixels in the upper part and 24x24 pixels in the bottom part. 
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M.G. Lofdahl: Evaluation of image-shift measurement algorithms for solar SH WFS 
Table 9: Results from seeing simulation. 



r (cm) 


LA 






No noise 










1% noise 






^ shirt V / 


cr (") 

*-* err V / 


(T lcT.\, 
w en / 17 tilt 


Fail (%) 




cr „ (") 

w shirt V. / 


cr (") 

err V / 


(T l(T,\, 
w efr/ ^ tilt 


Fail (%) 




5 


SDF 


l.n 


0.054 


0.049 


1.67 


1.011 


1.13 


0.158 


0.144 


1.22 


1.011 


5 


ADF 2 


l.n 


0.055 


0.050 


1.65 


1.010 


1.13 


0.169 


0.154 


1.27 


1.010 


5 


CFI 


1.12 


0.090 


0.082 


1.64 


1.012 


1.16 


0.181 


0.165 


3.44 


1.012 


5 


CFF 


0.84 


0.111 


0.101 


10.43 


0.832 


0.81 


0.509 


0.464 


2.06 


0.505 


7 


SDF 


0.84 


0.030 


0.037 


0.55 


1.009 


0.85 


0.115 


0.138 


0.52 


1.009 


7 


ADF 2 


0.84 


0.031 


0.037 


0.52 


1.009 


0.85 


0.123 


0.148 


0.54 


1.009 


7 


CFI 


0.84 


0.063 


0.076 


0.57 


1.010 


0.86 


0.132 


0.159 


1.27 


1.010 


7 


CFF 


0.69 


0.064 


0.077 


4.01 


0.859 


0.69 


0.235 


0.283 


5.57 


0.796 


10 


SDF 


0.62 


0.019 


0.031 


0.13 


1.009 


0.63 


0.095 


0.154 


0.25 


1.009 


10 


ADF 2 


0.62 


0.020 


0.032 


0.12 


1.009 


0.63 


0.102 


0.166 


0.28 


1.009 


10 


CFI 


0.62 


0.051 


0.083 


0.23 


1.009 


0.64 


0.110 


0.178 


0.54 


1.009 


10 


CFF 


0.53 


0.042 


0.069 


1.46 


0.871 


0.55 


0.179 


0.291 


2.51 


0.840 


1 c 
13 


SDr 


n a a 
U.44 


0.014 


U.U31 


U. 16 


1 nno 
1.UU8 


n a c 
U.43 


n noc 
U.U83 


n 1 nA 
0.194 


0.31 


i nno 
1.0U8 


15 


ADF 2 


0.44 


0.014 


0.033 


0.14 


1.008 


0.45 


0.091 


0.208 


0.33 


1.008 


15 


CFI 


0.45 


0.046 


0.104 


0.24 


1.008 


0.46 


0.099 


0.225 


0.48 


1.009 


15 


CFF 


0.38 


0.029 


0.066 


1.03 


0.876 


0.41 


0.157 


0.357 


1.06 


0.855 


20 


SDF 


0.35 


0.011 


0.033 


0.04 


1.008 


0.36 


0.082 


0.236 


0.20 


1.008 


20 


ADF 2 


0.35 


0.012 


0.035 


0.03 


1.008 


0.36 


0.088 


0.254 


0.22 


1.008 


20 


CFI 


0.35 


0.043 


0.126 


0.09 


1.008 


0.36 


0.095 


0.274 


0.37 


1.009 


20 


CFF 


0.30 


0.023 


0.067 


0.82 


0.878 


0.34 


0.149 


0.431 


0.72 


0.859 


5 


SDF 


1.11 


0.052 


0.047 


1.66 


1.010 


1.12 


0.111 


0.102 


0.97 


1.011 


5 


ADF 2 


1.11 


0.052 


0.048 


1.64 


1.010 


1.12 


0.119 


0.108 


0.93 


1.010 


5 


CFI 


1.11 


0.068 


0.062 


1.28 


1.012 


1.12 


0.121 


0.110 


1.03 


1.012 


5 


CFF 


1.00 


0.067 


0.061 


2.11 


0.919 


0.98 


0.193 


0.176 


3.10 


0.900 


7 


SDF 


0.84 


0.029 


0.035 


0.55 


1.009 


0.84 


0.078 


0.094 


0.39 


1.009 


7 


ADF 2 


0.84 


0.029 


0.035 


0.51 


1.009 


0.84 


0.083 


0.100 


0.38 


1.009 


7 


CFI 


0.84 


0.043 


0.052 


0.39 


1.010 


0.85 


0.085 


0.102 


0.39 


1.010 


7 


CFF 


0.77 


0.038 


0.046 


0.69 


0.930 


0.78 


0.137 


0.166 


0.87 


0.923 


10 


SDF 


0.62 


0.018 


0.029 


0.12 


1.008 


0.62 


0.063 


0.102 


0.11 


1.008 


10 


ADF 2 


0.62 


0.018 


0.030 


0.12 


1.008 


0.63 


0.068 


0.110 


0.13 


1.008 


10 


CFI 


0.62 


0.032 


0.053 


0.09 


1.009 


0.63 


0.070 


0.113 


0.10 


1.009 


10 


CFF 


0.58 


0.024 


0.040 


0.36 


0.935 


0.59 


0.114 


0.185 


0.34 


0.931 


13 


SDr 


n a a 
U.44 


0.012 


n noo 
U.U28 


U. 16 


1 nno 
1.UU8 


U.43 


U.U3D 


0.12/ 


n i n 
0.19 


i nno 
1.0U8 


15 


ADF 2 


0.44 


0.012 


0.028 


0.14 


1.008 


0.45 


0.060 


0.137 


0.19 


1.008 


15 


CFI 


0.44 


0.028 


0.063 


0.18 


1.009 


0.45 


0.062 


0.142 


0.18 


1.009 


15 


CFF 


0.41 


0.017 


0.038 


0.39 


0.938 


0.42 


0.102 


0.232 


0.36 


0.935 


20 


SDF 


0.35 


0.010 


0.029 


0.02 


1.008 


0.35 


0.053 


0.155 


0.07 


1.008 


20 


ADF 2 


0.35 


0.010 


0.029 


0.03 


1.008 


0.35 


0.057 


0.166 


0.09 


1.008 


20 


CFI 


0.35 


0.026 


0.075 


0.05 


1.009 


0.35 


0.060 


0.173 


0.06 


1.009 


20 


CFF 


0.32 


0.013 


0.039 


0.25 


0.939 


0.34 


0.098 


0.282 


0.27 


0.935 


5 


SDF 


1.11 


0.051 


0.047 


1.66 


1.010 


1.11 


0.083 


0.075 


1.08 


1.011 


5 


ADF 2 


1.11 


0.051 


0.047 


1.64 


1.010 


1.11 


0.087 


0.079 


1.05 


1.011 


5 


CFI 


1.11 


0.058 


0.053 


1.41 


1.011 


1.11 


0.087 


0.079 


1.06 


1.011 


5 


CFF 


1.06 


0.054 


0.049 


1.47 


0.965 


1.06 


0.129 


0.117 


0.88 


0.962 


7 


SDF 


0.84 


0.028 


0.034 


0.58 


1.009 


0.84 


0.055 


0.066 


0.38 


1.009 


7 


ADF 2 


0.84 


0.028 


0.034 


0.52 


1.009 


0.84 


0.058 


0.070 


0.34 


1.009 


7 


CFI 


0.84 


0.035 


0.042 


0.46 


1.009 


0.84 


0.059 


0.071 


0.36 


1.009 


7 


CFF 


0.81 


0.030 


0.036 


0.51 


0.971 


0.81 


0.091 


0.110 


0.41 


0.970 


10 


SDF 


0.62 


0.017 


0.028 


0.11 


1.008 


0.62 


0.043 


0.070 


0.05 


1.008 


10 


ADF 2 


0.62 


0.018 


0.028 


0.11 


1.008 


0.62 


0.046 


0.074 


0.06 


1.008 


10 


CFI 


0.62 


0.025 


0.040 


0.12 


1.008 


0.62 


0.047 


0.076 


0.05 


1.008 


10 


CFF 


0.60 


0.019 


0.031 


0.12 


0.973 


0.60 


0.075 


0.122 


0.15 


0.972 


15 


SDF 


0.44 


0.012 


0.027 


0.17 


1.008 


0.44 


0.037 


0.085 


0.16 


1.008 


15 


ADF 2 


0.44 


0.012 


0.027 


0.14 


1.008 


0.44 


0.040 


0.091 


0.13 


1.008 


15 


CFI 


0.44 


0.020 


0.045 


0.26 


1.007 


0.44 


0.041 


0.093 


0.15 


1.007 


15 


CFF 


0.43 


0.013 


0.029 


0.18 


0.974 


0.43 


0.067 


0.152 


0.23 


0.974 


20 


SDF 


0.35 


0.009 


0.027 


0.02 


1.008 


0.35 


0.035 


0.102 


0.03 


1.008 


20 


ADF 2 


0.35 


0.009 


0.027 


0.02 


1.008 


0.35 


0.038 


0.110 


0.03 


1.008 


20 


CFI 


0.35 


0.018 


0.053 


0.15 


1.007 


0.35 


0.039 


0.113 


0.03 


1.007 


20 


CFF 


0.34 


0.010 


0.029 


0.04 


0.974 


0.34 


0.064 


0.185 


0.11 


0.974 



Notes. Different CAs, IA always 2QI. The "Fail (%)" columns show the percentage of 4<x outliers, pi is the linear coefficient fitted to Eq. |7). The 
image scale is C.'41/pixel. In the top part, the FOV is 16x16 pixels or 6'.'6x6'.'6. In the center part, the FOV is 24x24 pixels or 978x9'.'8. In the 
bottom part, the FOV is 36x36 pixels or 14'.'8xl4'.'8. 
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M.G. Lofdahl: Evaluation of image-shift measurement algorithms for solar SH WFS 
Table 10: Results from seeing simulation. 



r (cm) 


LA 






No noise 










1% noise 






a mm (") 

^ shirt V / 


cr (") 

*-* err V ) 


(T lcT.\, 
w en / 17 tilt 


Fail (%) 




cr „ (") 

w shirt V. / 


cr (") 

err V / 


(T l(T,\, 
w efr/ ^ tilt 


Fail (%) 




5 


SDF 


1.09 


0.054 


0.049 


4.09 


1.011 


1.09 


0.114 


0.104 


3.50 


1.010 


5 


ADF 2 


1.09 


0.055 


0.050 


4.08 


1.011 


1.09 


0.122 


0.111 


3.47 


1.009 


5 


CFI 


1.09 


0.084 


0.077 


3.91 


1.011 


1.10 


0.135 


0.123 


3.88 


1.010 


5 


CFF 


0.88 


0.095 


0.087 


8.39 


0.859 


0.85 


0.212 


0.193 


8.50 


0.814 


7 


SDF 


0.84 


0.030 


0.036 


0.64 


1.009 


0.84 


0.079 


0.095 


0.53 


1.009 


7 


ADF 2 


0.84 


0.030 


0.036 


0.64 


1.009 


0.84 


0.085 


0.103 


0.53 


1.009 


7 


CFI 


0.84 


0.057 


0.069 


0.61 


1.010 


0.84 


0.096 


0.116 


0.57 


1.010 


7 


CFF 


0.71 


0.055 


0.066 


2.69 


0.879 


0.71 


0.135 


0.163 


3.20 


0.864 


10 


SDF 


0.62 


0.018 


0.030 


0.17 


1.009 


0.63 


0.064 


0.104 


0.21 


1.009 


10 


ADF 2 


0.62 


0.019 


0.030 


0.16 


1.009 


0.63 


0.069 


0.113 


0.22 


1.009 


10 


CFI 


0.62 


0.046 


0.075 


0.24 


1.009 


0.63 


0.080 


0.129 


0.19 


1.009 


10 


CFF 


0.54 


0.036 


0.058 


1.11 


0.888 


0.55 


0.107 


0.174 


0.82 


0.878 


1 c 
13 


SDr 


n a a 
U.44 


0.012 


n m"7 
0.02/ 


n no 


1 nno 
1.UU8 


n a c 


0.05 / 


n 1 on 
0.1 ly 


n i h 
0.1 / 


i nno 
1.008 


15 


ADF 2 


0.44 


0.012 


0.028 


0.07 


1.008 


0.45 


0.062 


0.140 


0.18 


1.008 


15 


CFI 


0.45 


0.041 


0.093 


0.12 


1.008 


0.45 


0.071 


0.163 


0.12 


1.008 


15 


CFF 


0.39 


0.024 


0.054 


0.89 


0.893 


0.40 


0.094 


0.213 


0.37 


0.884 


20 


SDF 


0.35 


0.009 


0.027 


0.06 


1.008 


0.35 


0.054 


0.156 


0.16 


1.008 


20 


ADF 2 


0.35 


0.010 


0.028 


0.05 


1.008 


0.35 


0.059 


0.169 


0.18 


1.008 


20 


CFI 


0.35 


0.039 


0.112 


0.07 


1.008 


0.36 


0.068 


0.197 


0.11 


1.008 


20 


CFF 


0.31 


0.018 


0.053 


0.88 


0.895 


0.32 


0.089 


0.256 


0.33 


0.886 


5 


SDF 


1.11 


0.051 


0.047 


1.50 


1.010 


1.12 


0.119 


0.109 


0.80 


1.010 


5 


ADF 2 


1.11 


0.052 


0.047 


1.48 


1.010 


1.12 


0.129 


0.118 


0.78 


1.010 


5 


CFI 


1.11 


0.064 


0.058 


1.16 


1.012 


1.12 


0.127 


0.115 


0.84 


1.012 


5 


CFF 


1.04 


0.059 


0.053 


1.33 


0.946 


1.04 


0.209 


0.190 


1.79 


0.933 


7 


SDF 


0.84 


0.029 


0.035 


0.40 


1.009 


0.84 


0.084 


0.102 


0.20 


1.009 


7 


ADF 2 


0.84 


0.029 


0.035 


0.39 


1.009 


0.84 


0.092 


0.111 


0.23 


1.009 


7 


CFI 


0.84 


0.040 


0.048 


0.30 


1.010 


0.84 


0.090 


0.109 


0.21 


1.010 


7 


CFF 


0.79 


0.034 


0.041 


0.34 


0.954 


0.80 


0.153 


0.184 


0.56 


0.950 


10 


SDF 


0.62 


0.019 


0.031 


0.08 


1.008 


0.62 


0.069 


0.112 


0.07 


1.008 


10 


ADF 2 


0.62 


0.019 


0.031 


0.08 


1.008 


0.63 


0.075 


0.122 


0.08 


1.008 


10 


CFI 


0.62 


0.030 


0.049 


0.24 


1.009 


0.63 


0.075 


0.121 


0.07 


1.009 


10 


CFF 


0.59 


0.023 


0.037 


0.11 


0.958 


0.60 


0.128 


0.208 


0.31 


0.955 


1 c 

1 J 




n a a 
U.44 


U.Ui4 


n mo 


n no 
U.UZ 


1 nno 




n nAO 
U.UoZ 


n 1 a 1 
U. 141 


n r\A 
U.U4 


i nn-7 
l.UU/ 


15 


ADF 2 


0.44 


0.014 


0.033 


0.02 


1.007 


0.45 


0.067 


0.153 


0.05 


1.007 


15 


CFI 


0.44 


0.026 


0.059 


0.38 


1.008 


0.45 


0.067 


0.152 


0.06 


1.008 


15 


CFF 


0.42 


0.017 


0.038 


0.04 


0.960 


0.44 


0.116 


0.263 


0.25 


0.958 


20 


SDF 


0.35 


0.012 


0.036 


0.01 


1.008 


0.35 


0.059 


0.171 


0.04 


1.008 


20 


ADF 2 


0.35 


0.013 


0.036 


0.01 


1.008 


0.35 


0.064 


0.186 


0.05 


1.008 


20 


CFI 


0.35 


0.025 


0.071 


0.44 


1.008 


0.35 


0.064 


0.186 


0.05 


1.009 


20 


CFF 


0.33 


0.015 


0.043 


0.02 


0.962 


0.35 


0.111 


0.321 


0.22 


0.959 



Notes. Different CAs, IA always 2QI. The "Fail (%)" columns show the percentage of 4<x outliers. p t is the linear coefficient fitted to Eq. {7}. 
The FOV is 24x24 pixels. In the top part, with an image scale of O'.'29/pixel, this corresponds to 6'.'9x6'.'9. In the bottom part the image scale is 
Or53/pixel and the FOV is 12'.'8x 12'.'8. 
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